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Abstract 

We investigate a double layer system with tight-binding hopping, 
intra-layer and inter-layer interactions, as well as a Josephson like cou- 
pling. We find that an antiferromagnetic spin polarization induces 
additional spin-triplet pairing (with S z = 0) to the singlet order pa- 
rameter. This causes an undamped collective mode in the supercon- 
ducting state below the particle-hole threshold, which is interpreted as 
a Goldstone excitation. 

PACS numbers: 71.45-d, 74.80.Dm, 74.50+r 

1 INTRODUCTION AND MODEL 

Collective density fluctuations in superconductors due to the breakdown 
of the global gauge invariance are well known theoretically Jl], ||. However, 
since these modes couple to charge oscillations, the long range Coulomb force 
usually pushes up their energies to the plasma frequency. One possibility to 
avoid the Coulomb interaction completely, is to consider spin fluctuations 
instead of charge fluctuations between the layers. In the following we will 
show the existence of such a sharp, collective spin mode in the gap, which 
might have been observedQ in inelastic neutron scattering on Y-Ba-Cu-O. 

We consider an electronic double-layer system described by the Hamil- 
tonian H = Hq + H$: 

H = Y, e k(4ka C lka + 4ka C 2ka)+ t ^4ka C lka + C \ka C 2ka) i 1 ) 
ka 
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HS — 2 Z^i ^ V\\ C ik+qa C ik'-qa' C ik'v' C ika ^ C ik+qa C jk'-qa' C jk'a' C ika 
kk'qcra' i 

+ J ( C lk+qa C lk'~qa' C jk'a-' C jka + C lk+qa C ]k'~qa' C ik'a' C jka)- (^) 

Here i& describes a tight-binding coupling between the two layers i = (1, 2), j = 
3 — i, while Vn (Vj_) are intra- (inter )-layer pairing interactions and the 
Josephson-like coupling J describes the coherent transfer of two particles 
from one layer to the other. In a previous publication ^ we treated this 
model using the Nambu formalism including vertex corrections to calculate 
charge fluctuations between the layers. 

In this paper we are primarily interested in the calculation of correlation 
functions involving the operator 

& = ^2 4fc| C 2fcT ~~ 4fcJ. C 2fc| ~ C lfcT C l/cT ~*~ C lfc| C lfc| 

k 

describing the difference of the spin polarization in the two layers and the 
operators coupling to it (A^. := c\ k ^c^_ ki ) 

$T = ~i Ek ^21 " A 12 " A 2i + A12, 

M = ~iEk 4fct C lfcT ~~ C lfet C 2fcT ~~ C 2kl C lkl + C \ki C 2k[- (^) 

The quantity M corresponds to the spin current between the two layers. 
<3?T and At describe pairing in different layers in a spin-triplet state with 
total spin S z = and are the real and imaginary part of the inter-layer 
triplet-pairing amplitude Aj_y := A12 — A21 ■ To shorten the notation, we 
introduce 

pij ■=J2 q 'k Dij ^k, D ij := *i® Tj , * k := (ci* T ,cI_ fci ,c 2fct ,4_ fc; ) t (5) 
k 

Ti ((Tj) being the Pauli matrices in the Nambu or two-layer space, respectively 
(examples: S = -P 30 , A T = -P 22 , $t = -P 21 )- 

2 ANALYTICAL RESULTS AND GOLDSTONE 
MODES 

In general the correlation functions <C pij^p lm ^> have to be determined 
numerically. However, for constant hopping tk = t with t,w < A (A is 
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the superconducting s-wave gap) and weak coupling the collective modes 
can be calculated analytically (for ujs^o <^ A) in the cases i (ii) of pure 
intra- (inter )-layer pairing. 



We obtain for case (i) case (ii) 

(It) 2 u? 



UJ 2 — LOg UJ 2 — COg ' 

iOn UJ V± — J 



«^ S>> = 2 ° 4 ^V-42A 2J 

w z — wfj A V|| + V± + 2 J 
«M,S» w -4ziV _ 2 -4iJVo-o r> 

W z — CJg W z — LOg 

uj 2 s = (2t) 2 + ujI (2t) 2 + cj ( 2 



! 

-{V\\ -V± + 2J) (2A) 2 -2J (2A) 2 



(Vl-J)(Vf,+J) JV Fi-J 2 iV ' 

(6) 

A spin polarization S with opposite sign in the two layers obviously induces 
inter-layer triplet-amplitudes At (3>t). 

These results are closely connected to the collective modes discovered 
in density-density-correlation functions like < P,P> (P := — P 33 ) in our 
previous work.Q For pure inter- (intra)-layer pairing one has the exact rela- 
tion 

<C P, P ^intor(intra) =< ^ & ^intra(inter) > (7) 

which follows from a unitary change of representation A = UAU>, \ ift) = 

U\ ip) with U := eaqp(-t7rEfc(cifci4fei + c 2fc|4fc|)) 

The Goldstone theorem || predicts the existence of excitations with 
vanishing energy, if a continuous, dynamical symmetry £1 ([Q, H] = 0) is 
spontaneously broken, e.g. the groundstate | 0) is not an eigenstate of Q. 
Thereby it can help to classify the resonances found in the correlation func- 
tions @ as so-called Goldstone modes connected with certain symmetries 
of H. 

Assuming pure singlet pairing in equilibrium, the superconducting ground- 
state is given by 

i =n( i +^n ei2eiiA iii, s +«^ 2e±A L5)) i o) (8) 
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Table 1: Broken symmetries if a» 7^ or a± 7^ 0. 



case 


parameters 


«|| / 


a x ^0 


1 


t, J, V\\ - v± = 


^03> ^33 1 ^13 > ^23 


^03, ^13; ^30> ^20 


2 


t,J^0,V\\ = v± 


^03 > ^13 


^03) ^13 


3 


t, J -» 0, Vf, 7^ VI 


^03, ^33 


^03 > ^30 


4 


*, J, V[| - VI ^ 


^03 


^03 



with Aug := An + A22 and A± t g := A12 + A21 being the singlet-order 
parameters for intra- and inter-layer pairing, respectively. The analytical 
calculations of case i (ii) refer to pure intra-(inter-)layer pairing with a± = 

(«|| = 0). 

Table|l| shows (for different parameters t, J, Vn , VjJ the symmetries fiy(<£) := 
exp(i(pP l:1 ), which are broken in the presence of intra- (inter )-layer pair- 
ing aii 7^ (aj_ 7^ 0). The spontaneous breakdown 1\)3 I 0\\,Q±) = 
\ 0\\ + 4>,9± + 4>) of the global gauge symmetry, which is generated by the to- 
tal particle number, in | 0ii , 9± } is a defining property of the superconducting 
phase (case 4), as it is invariably connected with non- vanishing Cooper-pair 
amplitudes ((Ay) 7^ 0). 

In Eq. H Goldstone modes (ug = 0) appear in case i (ii), if and only if 
t, J,V\\ — V± = (t, J = 0). We can identify these resonances in case i (ii) 
with the modes in table [l] in the cases 1 (3), which are connected with the 
symmetries ^23 (^30)- The transformations 

OmMIV-O = n fe ( 1 + Q ll e!2S " cos(2^)Aj s+axe^-LCA^g-s^^A^j,))^), (9) 
n 3O (<«|0||,0x) = n fe ( 1+Q! ll el2811 Aj |s +a ± e l2e ±(co S (20)A 1 LiS -ism(20)A 1 LiT ^|O) 

show that in both cases the S' 2 =0-component Aj_ t of the triplet-order pa- 
rameter is excited, which for O23 (^30) a l so creates non- vanishing expecta- 
tion values (At) ((^t)) and finite responses <C At, S S> (<C &t, S to an 
external spin polarization S. Thereby O23 mixes intra-layer pairs (Am g)with 
triplet-inter-layer pairs (A±^t), which is connected with a spin transfer be- 
tween the layers without breaking up Cooper pairs. On the other hand, 
Q30 leaves the modulus of the pairing amplitudes invariant, but creates a 
phase difference between A12 and A21, which is the origin of a supercurrent 
of inter-layer pairs, the so-called spin Josephson-effect. This terminology 
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is motivated by the close analogy to the usual Josephson effect, where a 
charge rather than a spin transfer is driven by a phase difference of intra- 
rather than inter-layer pairs. The density-modes O20 (^33)5 which corre- 
spond according to the relation (0) to the spin modes SI23 (^30) > can be 
observed as poles in <C P, P ^> calculated in our previous work Q rather 
than in <C S, S 3>. According to Eq. ^ all these modes cannot be excited in 
the absence of particle transfer (t, J = 0) between the layers. Finally, 
connects groundstates with different ratios of inter- and intra-layer-pairing, 
which for t, J = in case 1 and 2 are energetically degenerate. 

3 NUMERICAL RESULTS 



We carried out numerical calculations for < S,S > using two slightly dif- 



ferent effective masses 2mi/h 



1 eV _1 A- 



2m 2 /7V 



2 for 
0.25 



1.2 eV _1 A 

the two bands ± i& = ti z k' 2 /2m 2 /i and parameters: /i = 0.3 eV, uj c - 
eV (cut-off in k-space), (Vj| + V± + 2J)N = -0.44, (Vj| - V±)N = ±0.2 
(No: averaged density of states of the two bands). Then the pairing is mixed 
and for (Vn — V±)Nq < (> 0) dominated by intra- (inter)-layer pairing. 
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Figure 1: Im <C S, S ^> at T = in the case of dominant inter-layer (a) and 
intra-layer (b) pairing for different J. The spectral weight of the collective 
mode is given by the height of the <5-peak times 10 in (a) or times 100 for 
peaks marked with arrows and in (b). 



Fig. |H shows the imaginary part of the spin-polarization function for 
different J in the case a (b) of dominant inter- (intra)-layer pairing at T = 0. 
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The collective modes appear as 5-peaks below the onset of particle hole exci- 
tations around 60 meV. In the appropriate parameter range the resonances 
in case a (b) coincide with the poles calculated analytically in Eq. || in case 
ii (i), but they exist in a much larger parameter range. For larger positive or 
negative J-values than given in the figures the mode frequencies pass zero 
indicating an instability of the system. 

The spectral weight of the phase mode decreases for dominant inter-layer 
pairing going from negative to positive J as indicated by the spectral weight 
u>s in the approximation formula @. 

For positive J a further collective mode, the so-called amplitude mode, 
can be seen in both cases. It is inside the particle-hole spectrum for small J, 
but undamped for large J (arrows (a) or small peaks below the particle-hole 
threshold (b)). Because of the mixing of intra-layer and inter-layer pairing 
the spin excitation couples to both the phase <&t and the amplitude At of 
the triplet order-parameter. This causes two collective modes, which were 
already discussed in M. 
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Figure 2: Im <C S, S S> for different temperatures in the case of dominant 
inter-layer pairing for constant NqJ = 0.01 or constant NqJ = —0.3. 



The temperature dependence of the spin polarization function in Fig. 

2 shows the broadening of the collective mode at about 50 meV with 
increasing T. In case b a further collective mode appears at T c . 

To conclude, the anti-ferromagnetic spin polarization couples to the 
phase and amplitude of the triplet-order parameter with S z = 0. This 
causes a collective mode where spin-up and spin-down particle tunnel in 
opposite direction {spin Josephson-effect) . This might be connected with a 
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resonance found in magnetic neutron scattering on Y-Ba-Cu-0 at finite q, 
which shows the same temperature dependence as our mode[Q]. 
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